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We study the slippage on hierarchical fractal superhydrophobic surfaces, and find an unexpected 
rich behavior for hydrodynamic friction on these surfaces. We develop a scaling law approach for 
the effective slip length, which is validated by numerical resolution of the hydrodynamic equations. 
Our results demonstrate that slippage does strongly depend on the fractal dimension, and is found 
to be always smaller on fractal surfaces as compared to surfaces with regular patterns. This shows 
that in contrast to naive expectations, the value of effective contact angle is not sufficient to infer 
the amount of slippage on a fractal surface: depending on the underlying geometry of the roughness, 
strongly superhydrophobic surfaces may in some cases be fully inefficient in terms of drag reduction. 
Finally, our scaling analysis can be directly extended to the study of heat transfer at fractal surfaces, 
in order to estimate the Kapitsa surface resistance on patterned surfaces, as well as to the question 
of trapping of diffusing particles by patchy hierarchical surfaces, in the context of chemoreception. 

PACS numbers: 68.08.-p, 68.15,+c, 47.61.-k, 47.63.Gd 



I. INTRODUCTION 

Superhydrophobic (SH) surfaces have raised a consid- 
erable interest over the recent years [HQ ■ A SH surface is 
typically characterized by a contact angle exceeding 150°, 
which results from the combination of bare hydrophobic- 
ity and micro- and nano- scale roughness on the surfaces. 
This exceptional wetting property has motivated numer- 
ous studies to rationalize supcrhydrophobicity and device 
versatile SH materials, as well as to explore its implica- 
tion in the context of transport and fluid dynamics [H- 
[l2j . In particular, patterned SH surface were shown to 
exhibit low friction - superlubricating - properties in the 
Cassie state [l3|, i.e. when the liquid interface lies at 
the top of the roughness and thus strongly reduces the 
direct solid- liquid contact area [I3 - [l6j . This drag reduc- 
tion phenomena is associated with a large slippage of the 
fluid at the SH surface The quantity of interest 

is accordingly the effective slip length b g, entering the 
Navier boundary condition pj], [l8j for the averaged ve- 
locity field at the surface : b D sd z v = v w , where v w is the 
slip velocity at the wall, averaged over the lateral surface 

In this paper, we study the friction properties on mul- 
tiscale fractal SH surfaces. Fractal SH surfaces were in- 
deed shown to exalt super hydrophobicity [ljl [20| • Fur- 
thermore the fractal nature of the interface was shown 
to enhance the robustness of the Cassie state with re- 
spect to the Wenzel state (with the liquid impregnating 
the roughness), as compared to simple regular structures 
such as periodic stripes or pillars [2lU24|. However how 
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the fractal, multiscale nature of the interface impacts 
transport and hydrodynamic slippage remains largely un- 
known. Recently Feuillcbois et al. could obtain rigorous 
bound on the slip length on Hashin-Shtrikman fractals 
in a thin channel geometry [2o| . Another contribution 
showed that for a SH interface made of pillars, random- 
ness in the post structure does weaken the drag reduction 
as compared to periodic structure [26[. Also, recent ex- 
periments by Lee and Kim showed that a hierarchical 
surface with nanoposts lying on the top of microposts 
does not necessarily lead to an increase in slip length 
[27j . This points to the subtle connection between sur- 
face slippage and the structure of the underlying surface. 

A key difficulty is that there is not general analyti- 
cal theory describing hydrodynamic flows on composite 
surfaces. Therefore, analytical results for the slip length 
could only be obtained in specific and simple geometries 
with regular stripes [14l4l6|. while asymptotic calcula- 
tions on regular 2D post structures allowed to rationalize 
the limiting behavior of the slip length in the low surface 
fraction regime [28j. In this context, solving the Stokes 
equation with a fractal structure of boundary conditions 
appears hopeless. To bypass this difficulty, we resort to 
a scaling law approach for the fluid dynamics, which we 
put forward in a previous paper for regular surfaces [29| . 
In particular, this allows us to derive simple scaling laws 
to characterize the slippage properties on hierarchical - 
fractal - surfaces as a function of the roughness character- 
istics: surface solid fraction 4> s and fractal dimension df 
of the interface. The validity of this approach is then val- 
idated on the basis of numerical resolution of the Stokes 
equations. 

The main results of the present study are that: for 2D 
structures 

• The fractal dimension <if=l corresponds to a 
change of behavior of the friction: for fractal dimcn- 
lbocquet| s j on i ower than I , df < 1 , the effective slip length 
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diverges in the low solid fraction limit <f> s — > 0, while 
for df >1, the effective slip length saturates to a fi- 
nite value at low solid fraction. For df = 1 the 
divergence is logarithmic as 4> s — > 0. 

• For df < 1 the scaling approach predicts that the 
effective slip length behaves as 6 e ff ~ <t>^ a with an 
exponent a = (1 — df) / (2 — df) which depends di- 
rectly on the fractal dimension df of the interface. 
The explicit prefactor for this scaling behavior can 
be calculated analytically, leading to the prediction 
b c g ~ |i (^s) a for circular posts and valid in 
the low 4> s regime. This generalizes for fractal sur- 
faces the asymptotic result obtained by Davis and 
Lauga for circular posts on a square array (with 
a = 1/2) HI. 

• For a given solid fraction <f> s , regular surfaces (pe- 
riodic stripes or pillars) lead to higher slip lengths 
than for fractal surfaces. The difference between 
the two types of surfaces is enhanced at small solid 
fraction. 

• The value of the effective contact angle # e ff (charac- 
terizing the wetting properties) is not sufficient to 
quantify the amount of slippage on the surface: the 
effective contact angle is not a relevant parameter 
to infer the slippcrincss of a given surface. 

For ID structures, the slip length for hierarchical frac- 
tal surfaces does saturate to a finite value at small surface 
fraction <p s and is therefore always far below the result 
for regular stripes. 

Here we consider hierarchical fractal SH surfaces ob- 
tained by the periodic repetition of a unit cell. We denote 
by L the period of the structure. We make an iterative 
geometric construction of the unit cell. For the first it- 
eration (n=l) the unit cell is divided in p parts in each 
direction among which k patches represent a solid sur- 
face, the other parts corresponding to the gas interface 
with negligeablc friction. At all further iterations n > 1, 
each solid patch is replaced by the pattern defined at the 
first iteration, with the initial period L replaced by the 
smaller scale of the solid patch. On Fig Q] we have repre- 
sented two examples of the unit cells of hierarchical SH 
surfaces that we have studied in 2D (a) and in ID (b) at 
the first iteration n=l (al and bl) and at the second iter- 
ation n=2 (a2 and b2). Using this recursive construction 
for the hierarchical multi-scale structure, there is self sim- 
ilarity only over a finite range of scales, the "true" fully 
fractal surface being obtained in the limit n — > oo. The 
fractal dimension for such structures is: 



d f = 



log k 
\ogp' 



(1) 



The initial solid fraction (for n=l) <f>o, is equal in fD 



to 0o = -, and in 2D to <po — ^ (here patches with a 



squared geometry are considered). The solid fraction at 
the second iteration ?i=2 is then: (f> 8 2 = (</>o) 2 and more 
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(32) 



□ 



(bl) 



ID 

k=2, p= 
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FIG. 1. Iterative construction of the hierarchical surfaces 
considered in this study, in 2D (al, a2) and ID (bl, b2). 
Filled patches correspond to the solid-liquid interface, while 
open patches correspond to the gas-liquid interface. The lo- 
cal boundary condition are no-slip on the filled patches and 
perfect slip on the open patches. In 2D for k=5, p=3, (al) 
first iteration n=l and (a2) second iteration n=2. In ID for 
k=2, p=3: (bl) n=l and (b2) n=2. 



generally the solid fraction at an iteration n is: 



II. A SCALING LAW APPROACH 

We first present the basis of the scaling law approach 
which we developed to evaluate the slip length on SH 
fractal surfaces. The starting point is a hierarchical re- 
lationship which we proposed in Ref. f29j . and showing 
that for a SH surface of roughness periodicity L, a (small) 
solid fraction <j) s and slip length b s on the solid at the top 
of the roughness, the effective slip length b e s takes the 
following expression: 



b c s ^ b Q + — . 



(2) 



In this equation, bo stands for the slip length for the same 
geometry when there is no slip on the solid at the top of 
the roughness (b s = 0) f3C}. Typically we showed that 
bo r~j a/<p s where a the typical lateral size of the roughness 
feature. 

We briefly recall the physical idea underlying this rela- 
tionship. It is obtained by identifying the averaged fric- 
tion force on the surface, which is by definition F s = 
ArjU/bes - with A the lateral area and U the slip veloc- 
ity -, with the total viscous force acting on the top of 
the solid roughness (say, posts or stripes), which writes 
F n = r)<f> s (7), with (7) the average shear-rate on a sin- 
gle roughness feature. As proposed in Ref. [2!|, a sim- 
ple estimate for this quantity is (7) ~ U/(a + b s ), with 
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a the typical size of the solid roughness feature. This 
force balance leads to b e g ~ a/4> s +b s /<f) s , and thus to 
Eq. ([2]) . Note that this derivation neglects hydrodynamic 
interactions between different roughness features and is 
expected to be valid in the limit of small solid fraction 
4> s . In Ref . |29[ , we have checked in detail the validity of 
the above analysis and of Eq.®. 

Here we extend this scaling analysis to the case of hi- 
erarchical fractal SH surfaces: the main idea is to use 
Eq.@ as a hierarchical relation between the different 
scales of the roughness, with b s characterizing the slip- 
page of the patch at lower scales. 

A. Wettability and contact angle 

Let us first consider the wettability of these surfaces. 
This is usually characterized by the effective contact an- 
gle # e ff of a droplet on the surface. Using Cassie-law, 9 e s 
is obtained in terms of ratio between the solid-liquid and 
gas-liquid contact area. Introducing <f> s the solid surface 
fraction, one has 

cos# c ff = <p s (l + cos6* ) — 1 (3) 

with 9q the contact angle on the bare, smooth surface 
with the same chemical caracterisitics. As pointed above, 
the solid fraction is given by 4> s = 4>si = <fio a t step 
n = 1, while at step n = 2, (ft s = <p S 2 = 4%. For a 
surface characterized by a hierarchy of n successive scales, 
then (j) s = <j) sn = 0q. Going up in the hierarchy, the 
solid fraction (j> s thus decreases and the contact angle on 
those hierarchical fractal surfaces converges exponentially 
towards 180°: 

cos6» e ff + l -> (4) 

n— >oo 

Let us now turn to the slippage properties. 

B. Slippage on a two scale roughness 

To illustrate the approach, we begin by studying the 
simplest hierarchical surface obtained with two steps in 
the iteration process (n=2), i.e., two roughness scales 
as represented in FigfTJ (a2) and (b2). Such hierarchical 
surfaces are now accessible from experiments thanks to 
the progress in nanofabrication techniques [2l|, [23| . 

Let us introduce b\, the effective slip length obtained 
after one iteration n=l, sec Fig[JJ (al) and (bl). A key 
remark is that, when going from the first n — 1 to the 
second step n~2, the higher-level structure inside each 
main solid patch is equivalent to the initial bare structure 
at n = 1, up to a rescaling of the size by 1/p (p being the 
size ratio between two scales, as defined above). 

Using Eq[51 we then deduce that the effective slip 
length is equal to: 

b 2 = b 1 + (5) 



where we recall that tfio is the solid fraction at the ini- 
tial scale, in FigUJal)-(bl). For ID structures like in 
FigOJbl), 0o = k/p, so that one obtains 

&a = &ix(l + i), (6) 

For 2D structures like in FigUJal), <fio = k/p 2 and 

&2 = 6ix(l + |), (7) 

The slip length is obviously higher at the second itera- 
tion. There is an increase in slip length for hierarchical 
"fractal" SH surfaces and this increase is larger for 2D 
structures than ID. Such a difference between the ID 
and 2D structures was already pointed out for regular 
periodic SH surfaces (29j . 

C. Slippage on hierarchical fractal surfaces 

Now, "real" SH fractal surfaces are composed by an 
infinite number of iterations, or at least by a large number 
n of it, so that the fractal dimension is defined over a 
broad range of scales. To estimate the slip length, we 
extend the previous analysis for any number of iterations. 

At a subsequent iteration n + 1, the solid patches are 
divided once again into p parts, see FigfT] As an alter- 
native view, one may consider that the system at the 
subsequent iteration n + 1 is made by first assembling k 
systems obtained at the previous iteration n in order to 
gather a structure with the correct geometry, and then 
perform a global rescaling by 1/p to obtain the desired 
structure at step n+1. Accordingly, the structure inside 
the largest patches at the iteration n+1 can be seen as 
the previous structure at iteration n, up to a rescaling 
factor p between scales. 

One may deduce directly the recursion relationship for 
the slip length from the above procedure. At scale n+1 
the slip length over each of the largest patches is then 
b n /p (since we have gathered structures from iteration 
n and rescaled by 1/p; remind that the slip length does 
scale with the global lateral size of the system). Now 
using Eq. [5J the slip length 6 n +i at iteration n + 1 can 
be deduced as a function of b n as: 

b n+ i = 61 + ^ x ^ (8) 
<Po P 

The last term thus gathers all the slippage effects on all 
above scales > n in the hierarchy; 61 is the slip length of 
the elementary structure, as in FiglTJal)-(bl). This can 
be rewritten as 

bn+i = bi + b n x with x = — — . (9) 

P<Po 

so that 

n — 1 n 

b n = hx y^ = 5!— — (10) 

f— i 1 — X 
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Depending on the parameter x, different scenarios oc- 
cur: 



• For x < 1, the sum in Eq.(|10p converges exponcn- 
tiaily to a finite value in the limit n — > oo: 



61 x 



1 



1-x 



(11) 



• For x = 1, the effective slip length at an iteration 
n is proportional to n, b n = n X 61. Now, replacing 
the scale index n by the corresponding solid surface 
4> s at the corresponding scale, S = 0q so that n = 
log 0s /log (^0, then one gets 



&eff(0s) = = h X 



l0g( 



log 0o ' 



(12) 



The slip length thus diverges logarithmically with 
the solid fraction 0, in the limit of small 0,. 



For x > 1, the largest term in Eq. (jT0| dominates 
the effective slip length, so that the slip length 
has the asymptotic expression in the limit n — > 00 
(0 S 0) 



Ocffl 



(13) 



with n 



log 00 ' 



As a side note, one may remark that these results could 
be generalized to the situation with a finite slip length 
60 on the liquid-gas interface. Our present results thus 
propose an upper bound for the effective slip length in the 
latter situation. We leave this generalization for future 
work. 



maximal amplification by a factor of two for the slip 
length as compared to a regular SH surface defined 
with only one roughness scale. In this regime, the 
slip length does not diverge in the limit of perfect 
superhydrophobicity, 6 e g = 180°, or S = 0. 



• The peculiar case df=0 (x = 1, i.e. k=l) corresponds 
to a flow parallel to stripes. We recover in this case the 
slow logarithmic divergence of the slip length with the 
solid fraction b n 
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log 0J as obtained in Refs. 
This result is non trivial as in the present approach we 
do not cxpliclity take into account the hydrodynamics 
interactions between stripes. 



2. 2D structures 

In 2D, a variety of situations occurs depending on the 
ration x = p/k. The latter parameter quantifies the 
dense or sparse character of the solid filling on the sur- 
face. 

• For 1 < df < 2 (x = p/k < 1), the solid filling is quite 
dense - this is the situation depicted in FigJIJal)-(a2) 
-. In this regime, the slip length of the SH hierarchical 
surface saturates to a value 
low S . 

• For df = 1 (x = p/k = 1), we find a logarithmic 
divergence of the slip length with the solid fraction b n = 
n x bi = b\. We find the same divergence as for a 



foifeZ^ m the limit of 



log<p 

ID regular SH surface of stripes. 

• For df <1 (x — p/k >1), the solid filling is sparse. In 
this case the effective slip length diverges as a power law 
of the solid fraction. Using Eq[T31 we get for large n, i.e. 
small S : 



D. Discussion: 1D/2D hierarchical surfaces versus 
regular surfaces 

We now summarize the different scaling expression for 
the slip length for ID and 2D structures and for the dif- 
ferent values of x. A key quantity is the fractal dimension 
of the structure, df = . 

• For ID structure, x = 1/k and k < p, so that x < 1 
and df <1. 

• For 2D structure, x = p/k and k < p 2 so that one finds 
three different cases: for x <1, 1 < df < 2; for x = 1, 
df=l; and for x >1, df <1. 



ID structures 



1 



"effV 



x- 1 



bi_ x S 



(14) 



with again df — j^S— the fractal dimension. 

We can compare the scaling behavior predicted in 
this last regime with the one obtained for a SH surface 
made out of a bidimensional regular pattern of posts for 
which the slip length behaves for small S as b c ^ ~ 

< i, this shows that for a given solid 



Since 



1- df 

2- df 



fraction S , the slip length on the periodic patterns of 
post is larger than of 2D SH fractal surfaces. This is 
a counter- intuitive result, which shows that fractality, 
though strongly increasing the non- wettability, is far less 
efficient in enhancing dynamic properties. 



In ID, two situations occur: 

• For 0< df <1 (x <1), then = friary, and the 
maximum value of b x /b\ is 2 (for k = 2). This shows 
that the hierarchical fractal geometry only generates a 



E. Alternative scaling law approach in the dilute 
regime 

Here we propose a simpler, alternative description to 
account for the scaling laws above. We focus on the case 
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of a sparse solid filling in 2D, corresponding to dt < 1 
and x larger than unity. 

Then in the dilute regime, when the hydrodynamics 
interactions can be neglected, the effective slip length 
6 e ff can be estimated directly from the calculation of the 
viscous friction force F of a liquid of viscosity 77 on a unit 
cell of size L of the hierarchical SH surface at a given 
step n. We consider a flow leading to a slip velocity U at 
the surface. By definition of the slip length, the friction 
force on the surface takes the expression: 



F = T) 



U_ 

b c ff 



L 



(15) 



with T] the shear viscosity of the fluid. Now the friction 
force can be calculated as the sum of the forces on each 
solid post. Separating the scale of the post from the 
lateral scale between the post, this friction force can be 
calculated from the Stokes equation applied to the posts 
considered as individual stationary solid object in a fluid 
moving at velocity U far from it: 



F v = ^2 XV a nU 



(16) 



solid patches 



Here the sum runs over all the solid patches on the sur- 
face, with characteristic lateral size a n at the smaller 
scale n; here x is the numerical prcfactor for the Stokes 
law. The latter depends only on the geometrical charac- 
teristics of the solid post. For example for a solid disk 
with radius a, one has x — \ x if El- The factor 4^ 
comes from the Stokes flow calculation around a disk, as 
calculated in Ref. [3l[ , while the 1 /2 prefactor stems from 
the fact that only the upper-half of the disk undergoes 
friction. 

Now at an iteration n, the lateral size of the solid post 
scales like a n ~ L/p n and there are k n of such posts. 
Altogether equating the two expressions for the friction 
force, we thus obtain the scaling law for the slip length 
in the dilute regime: 



Soft 



(17) 



Using the definition of i 
behavior 



and df this leads to the scaling 



'off 



(18) 



and we recover the same result as in Eq ll4[ obtained with 
the previous hierarchical approach. 

Now we finally remark that further insight can be ob- 
tained from the above analysis to predict the prefactor 
of the previous scaling behavior. We will focus on posts 
with disk shape for which the prefactor x m the viscous 
force, in Eq. (fT6"|) . can be exactly calculated, as pointed 
out above. 

We first start by noting that applying a similar analysis 
to a regular square pattern of disk-like posts yields b e s = 
3/16 x L 2 /a with a the disk radius, and we used x = 16/3 



for the disks as discussed above. Using <f> s = ira 2 / L 2 for 
the surface fraction, one gets the following expression for 
the slip length 



freff 



3^ L 
"16" 



(19) 



The prefactor ^5 is the exact one calculated by David 
and Lauga [H] using asymptotic expansions and con- 
firmed in numerical resolution of the hydrodynamic prob- 
lem [2£| . This demonstrates the power of the scaling 
approach in deriving relevant expression for slip lengths. 
The above analysis can also immediately extended to pre- 
dict the slip length on an array of posts of any shape (e.g. 
ellipsoids) and any pattern symctry. 

Coming back to the hierarchical pattern, a similar 
analysis for disk- like posts allows to show that b e g/L = 
2/x x x n with x = 16/3 for disks. Using S = \ (k/p 2 ) n , 
where the 7r/4 factor takes into acount the disk shape of 
the posts, this equation can be rewritten as 



b c s : 
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-L 



(20) 



The previous result for a regular square pattern of disk 
posts is exactly recovered by putting df = in this equa- 
tion. 

Furthermore this result can be rewritten as 6 c ft = 
b\ x a;"" 1 , with b\ = ^x: this expression is identical 

to Eq. (|T3"1) obtained by the alternative iterative deriva- 
tion, for x 1 (sparse structure). This demonstrates 
the equivalence between the direct and hierarchical cal- 
culations for the slip length on fractal surfaces. 



III. NUMERICAL RESULTS 

In order to assess the validity of the above scaling ap- 
proach, we perform a numerical resolution of the full 
hydrodynamic equations complemented by the compos- 
ite boundary conditions on the hierarchical surface, as 
sketched for example in FiglTJ Specifically, the hydrody- 
namic velocity field obeys Stokes equation in the bulk, 
with a no-slip boundary condition on the solid patches 
and a perfect slip boundary condition on the remaining 
liquid-gas interface (open patches in FigO}. 

To this end, we follow the numerical approach devcl- 
opped by Cottin-Bizonnc et al. in Ref. [32j . We only 
recall here the basic steps of this approach, and a more 
detailed description can be found in [32j. A shear flow 
is considered over a composite surface characterized by a 
heterogeneous slip length pattern. The boundary is mod- 
eled by a pattern of local slip lengths on a planar surface. 
The hydrodynamic equations for flow profile are rewrit- 
ten analytically in terms of a boundary integral method, 
which is then solved numerically. The effective slip length 
is deduced from the asymptotic flow profile far from the 
surface. 
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Using this numerical approach, we calculate the effec- 
tive slip length b e g for SH hierarchical fractal surfaces 
as the ones represented on FigQ] The computed slip 
lengths b c s/L will be calculated for different SH hierar- 
chical surfaces, for both ID and 2D structures, as well as 
for different values of k and p and iteration steps, n. The 
numerical values will then be compared with the predic- 
tions obtained using the scaling law approach. 

In practice, due to the large range of spatial length 
scales involved, the numerical calculations require a 
higher discretization than the one performed with a 
regular square lattice of solid patches as considered in 
[29|. The calculations involve subsantially larger matri- 
ces, which limits the range of fractal surfaces that can 
be considered. Consequently, the numerical method is 
limited to only a few iterations as shown in FigJ3] and 
[5l This will however prove sufficient to explore the main 
tendencies predicted by the hierarchical approach. 



Numerical results for slippage over hierarchical 
2D structures 



We have first plotted in Fig f2] the numerical results for 
the effective slip length b a s/L as a function of the solid 
fraction <j> s , for SH surfaces made out of both periodic 
posts on a square lattice and hierarchical fractal surfaces 
obtained with different parameters k and p. A first strik- 
ing result emerging from this plot is that the surfaces 
with regular structures of posts are far better than all 
the fractal surfaces in terms of slip length for a given 
solid fraction. 

The detailed comparison between the numerical results 
and the incremental scaling laws is then shown in Fig|3] 
A general conclusion from this plot is that the scaling 
predictions in the previous section are in good agreement 
with numerical results. 

(i) Case x — p/k < 1 i.e. fractal dimension df >1: 
two cases have been solved numerically with hier- 
archical structures corresponding to {fc,p}={5,3} and 
{k,p}={6,3}. At the first iteration n = 1, these struc- 
tures have a surface fraction of <pQ = 0.55 and <po = 0.66 
respectively. For further iterations, the solid fraction 
obeys </> s (ra) = 

The scaling predictions in this case suggest a finite slip 
length, which saturates at small cj> s . This is indeed the 
tendency shown by the numerical results. 

To be more quantitative, we use the results for the 
slip length in EqUHl b n = h(l - x n )/(l - x) with n = 
log(l/0 s )/log(p 2 /fc). One has to fix the numerical pref- 
actor b\ in the previous prediction: here we choose to 
calculate this value numerically for a surface structure at 
step n = 1, such as on FigQTal). Using this calculated 
parameter as an input, a quantitative agreement is then 
found between the scaling prediction and the numerical 
results, as shown in Figf3] 

(ii) Case x = k/p = 1, i.e. fractal dimension df =1: the 




FIG. 2. Numerical results for the normalized effective slip 
length fteff /ion 2D structures, as a function of the solid frac- 
tion <f> B (in log scale) for posts (+) and hierarchical surfaces 
(♦:{fc,p}={2,3}, ■= {fc,p}={3,3}, o:{fc,p}={5,3}). The dot- 
ted line corresponds to the scaling law for regular posts from 
[2j|. Typical error bars are represented for {fe,p}={2,3}. 



numerical calculations have been performed for a hierar- 
chical structure with {fc,p}={3,3}. At the first iteration 
n = l, this structure has a surface fraction of 4>o = 0.33. 
Note that only four iterations could be calculated numer- 
ically. 

The scaling laws predict in this case a logarithmic 
dependence on the solid fraction <p s and we thus com- 
pare the numerical expression with the estimate b n ~ 

logfo- 1 )- 

As above, we compute the numerical prefactor for this 
behavior by calculating the value of the slip length for 
a given iteration. We calculated numerically b\ at the 
7i = 1 iteration (with <j> s i = 0.33) to fix the prefactor of 
this <j) s dependence. The predicted slip length takes the 
following expression in this case b c g(4> s ) = b\ X n, with 

71 = log(0 s )/log(0o). 

Using this value as an input, we obtain a quantitative 
agreement between the predictions and the numerical re- 
sults for all (f> s < 0.1, demonstrating the robustness of 
the scaling approach up to relatively high solid fraction. 

(Hi) Case x = p/k >1, i.e. fractal dimension df <1: 
the numerical calculations have been performed for a 
hierarchical structure with {fc,p}={2,3}, with df = 
logfc/logp = 0.63. At the first iteration n = 1, this 
structure has a surface fraction of 4>o = 0.22. Note that 
only three iterations could be calculated numerically for 
this very sparse system. 

In the low <f> s regime, the hierarchical approach pre- 
dicts a slow divergence vs surface fraction, b Q g ~ <fi~ a 
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FIG. 3. Comparison between scaling law predictions and 
numerical results for 2D structures: normalized slip length 
bes/L as a function of the solid fraction cf> s (in a log-log plot) 
for hierarchical fractal surfaces with different fractal dimen- 
sions df. Theoretical predictions of the scaling approach, 
Eq. (|10[) . are shown as lines for the various fractal dimen- 
sions df under consideration (see text for details). Numeri- 
cal results are from top to bottom: (i) df <1: {k,p}—{2,3} 
(♦); (ii) d f =1: ■{fc,p}={3,3}; (iii) d f >1: o{fc,p}={5,3} 
and A{fc,p}={6,3}. The top dotted line is the predicted 
asymptotic scaling in Eq. (|20[) , b c g ~|L(^0 S ) a with a = 
(1 — df)/(2 — df). Here a ~ 0.27 for the fractal dimen- 
sion df — 0.63 corresponding to the hierarchical surface with 
{fc,p}={2,3}. 



eter, the agreement is slightly less good, as the scaling 
approach gives only a fair, semi-quantitative estimate to 
the slip length b\. 

Altogether the numerical results do confirm the theo- 
retical predictions. A quantitative agreement can even 
be obtained by fixing the prefactor of the scaling behav- 
ior, using the numerical value of the slip length of the 
elementary structure (at the iteration n — 1 or n = 2). 

Finally we consider in this paragraph the relationship 
between slippage and contact angle. The hierarchical 
fractal SH surface is characterized by its effective contact 
angle as defined in Eqf5] for a given iteration number n. 
We have plotted in Fig|4] the effective slip length as a 
function of 9 c s for both hierarchical fractal SH surfaces 
and for SH surfaces made out of periodic posts (numerical 
results and scaling laws). Here we defined the effective 
contact angle as # e ff = acos(0 s — 1), thus assuming for 
illustration purposes that the contact angle on the corre- 
sponding smooth surface is equal to 9q = 90°. 

This plot illustrates immediatly that the slip length is 
not a unique function of the contact angle 9 c s: in con- 
trast, it depends strongly on the geometry of the under- 
lying structure for a given value of 9 c g , as characterized 
here by various fractal dimensions d / . The hydrodynamic 
properties of a SH surface and thus of a hierarchical frac- 
tal SH surface can not be simply determined from its 
effective contact angle. A strongly non-wettable surface 
with a contact angle of 9 e g ~ 180° may well exhibit a 
low slip length and thus be inefficient in terms of drag 
reduction. 



with a = (1 — df)/(2 — df) ~ 0.27. To be more precise, 
we show in Fig[3] (top dotted line) the analytical asymp- 
totic prediction obtained in Eq. (|20"]) . 5 c ff ~ |L (j<^s) " 
with a = (1 — df)/(2 — df). This expression provides a 
result in good agreement with numerical results for small 
4>s [Hjj. Note however that in this case, the asymptotic 
(algebraic) scaling regime is still not fully reached for the 
smallest solid fraction considered. As can be seen from 
Eq. (fl"U| , the speed of convergence is actually fixed by the 
parameter x = p/k, which accounts for the sparsity of 
the structure: convergence to the scaling regime is ac- 
cordingly faster for more sparse structure with x 3> 1, 
i.e. small df, while x = 3/2 in the present case. 

However, as we now show, a quantitative agreement 
can even be reached using the full expression in Eq llOl to 
evaluate the slip length. We use the numerical value for 
62 at the n = 2 iteration (with (f> S 2 — 0.05) as an input 
parameter to fix the prefactor for this law in order to 
compare to numerical results. With this reference value, 
the slip length takes the following expression 6 c ff(0s) = 
6 2 x (x n - l)/(x 2 - 1), with n = log(0 s )/log(0 o ). As 
shown in Fig[31 this prediction is in quantitative agree- 
ment with the numerical results. Note that if one uses the 
n = 1 iteration as a reference for the prefactor param- 



B. Numerical results for slippage over hierarchical 
ID structures 

We finally turn to surfaces with ID structures. We 
have plotted on Figf5]thc numerical results and the scal- 
ing laws for b e f{/L as a function of <p s (in log scale) for 
ID hierarchical fractal SH surfaces characterized by dif- 
ferent couples {k,p}. We have also plotted the analytical 
expression expected for a flow parallel to periodic grooves 
[14J. As for the 2D structures, the fractal surfaces yield 
values for the slip length which arc far smaller than the 
one obtained on regular grooves. Furthermore as seen on 
this figure, we find a good agreement between the numer- 
ical results and the results of the scaling law approach. 
As for the 2D structure above, the prefactors for the scal- 
ing relationships are obtained by matching the slip length 
for the n = 2 or n = 3 iteration, as the scaling approach 
is expected to be valid for relatively small values of 4> s . 

IV. CONCLUSION 

In conclusion, we developed scaling laws for slippage 
on hierarchical fractal SH surfaces, providing explicit ex- 
pressions for the slip length as a function of the solid 
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FIG. 4. Normalized effective slip length b e g/L in logarithmic 
scales, as a function of the effective contact angle # e ff for dif- 
ferent 2D SH surfaces. Symbols are the same as in Fig(3] The 
black dotted line (top) is the scaling law for SH surface made 
out of periodic posts as predicted from Ref. [29^ . 
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FIG. 5. Comparison between the numerical results and the 
scaling law approach for ID structures: numerical results for 
the normalized slip length b a s/L versus solid fraction <f> a . 
Theoretical predictions of the scaling approach, Eq. (|10|l . are 
shown as lines for the two systems under consideration (see 
text for details). Numerical results are from top to bottom: 
{fc,p}={2,4} (■); {fc,p}={2,3} (♦J.The dash-dot line is the 
analytical solution for a flow parallel to periodic grooves as 
predicted in Ref. [l4T |. 



fraction <j> a and fractal dimension df of the surface. For 
two dimensional structures, the slip length is found to 
scale as 6 c ff ~ 4*7° with a = |— r an d df the fractal 
dimension (here df = log k/ logp for the hierarchical sur- 
faces under consideration). 

We checked the validity of this approach by numerical 
calculation of the full hydrodynamic equations. The nu- 
merical results show that the simple scaling description 
remains valid up to relatively high values of solid frac- 
tions, even though the scaling laws do not take properly 
into account the hydrodynamic interactions. A quanti- 
tative agreement can be obtained by fixing the prefactor 
of the scaling behavior, using the numerical value of the 
slip length of the elementary structure (at an iteration 
n = 1 or n = 2). For df < 1 the scaling laws provide 
a good prediction for the slip length up to 4> s = 0.1 and 
for df > 1 for which there is a fast saturation of the 
slip length, the agreement is good even for higher values 
of (j) a . Globally the numerical results confirm that the 
scaling approach captures the physical mechanisms at 
play: although hydrodynamic interaction are neglected 
between posts, scaling results are shown to be predictive 
over a broad range of solid surface fraction. 

Altogether, our study shows that for a given solid frac- 
tion (f> s , fractal SH surfaces are less efficient than SH sur- 
faces with regular pattern for drag reduction purposes. 
Furthermore for a given effective contact angle # e ff, the 
slip length on fractal SH surfaces decreases with increas- 
ing fractal dimension. This shows that slippage depends 
strongly on the geometry of the underlying structure for 
a given value of # e ff and the amount of slippage on a 
fractal SH surface can not just be determined from its 
effective contact angle. A strongly non-wettable surface 
with a contact angle close to 180° may well exhibit a 
low slip length and thus be totally inefficient in terms of 
drag reduction. This result is actually counter-intuitive: 
it shows that a surface with vanishingly small solid frac- 
tion leads to a finite friction force. 

The use of fractal SH surfaces appears however as a 
good compromise between the positive effect of stabiliza- 
tion of the Cassie state and a moderate reduction of the 
slippage effect. Very large slip length and corresponding 
large drag reduction can indeed be observed on regular 
SH surfaces, but usually at the expense of a low stabil- 
ity of the Cassie state versus pressure. A recent study 
has demonstrated the recovery of the Cassie state after a 
transition towards the Wenzel state [34| ■ It would be thus 
interesting to consider surfaces with fractal structure in 
the context of stabilization issues [22[ . 

Finally we mention that all results found in the present 
study for slip length can be immediately extended to 
heat transfer across fractal surfaces. In this case, the 
role of 6 G ff is played by the Kapitsa resistance [H, HU , 
which quantifies the temperature slip at the surface be- 
tween two different materials. A key point indeed is that 
the scaling arguments proposed here rely ultimately on 
the Laplacian nature of Stokes equation, so that all scal- 
ing results presented here do apply equally to any trans- 
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port property obeying a Laplacian equation, such as heat 
transport. Accordingly, one expects that the Kapitsa re- 
sistance (and corresponding Kapitsa length) at the inter- 
face with a hierarchical, fractal surface will obey similar 
scaling laws versus the solid fraction as the ones found 
here for the slip length: the Kapitsa length is expected to 
scale as 4>J a with a = [1 — df) J {2 — df) for surfaces with 
fractal dimension df < 1 in the low (f> s regime, while 
it is expected to saturate at a finite value for df > 1. 
This scaling results suggest therefore that for fractal di- 
mensions df > 1, it may be possible to have a relatively 
low Kapitsa surface resistance, thus maintaining a good 
heat transfer between the solid surface and the liquid, 
while keeping a strongly non-wetting superhydrophobic 
behavior. Such a compromise would be valuable for heat 
transfer problems in cooling pipes. 



As an extension of this study, our results could also be 
applied in a different field to the question of trapping of 
diffusing particles by patchy surfaces |36j . which is rele- 
vant to the questions of chemoreception and ligand bind- 
ing to surface receptors [33, Hll . A counter- intuitive pre- 
diction emerging from our analysis is that in the regime 
of an effective fractal dimension f < d c g < 2, a hierar- 
chical surface with vanishingly small surface fraction of 
receptors may yet result in a finite effective trapping rate. 
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